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SUMMARY 


This  report  provides  a demonstration  of  an  interpolative  model  for  crack  propagation  life 
predictions  at  elevated  temperature.  The  model  is  based  on  the  hyperbolic  sine.  Two  spiecimen 
geometries  were  tested  under  stress-temperature-time  spectra  that  represent  advanced  gas 
turbine  disk  operating  conditions.  A critique  of  model  predictive  ability  is  presented. 
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SECTION  I 
INTRODUCTION 

I 

I 

Historically,  the  methods  for  predicting  low-cycle  fatigue  (LCF)  life  have  produced 
conservative  underestimfetions  of  total  useful  life,  resulting  in  costly  early  replacement  of  LCF 
limited  gas  turbine  engine  rotating  components.  Accurate  total  LCF  life  predictions  must 
consider  (1)  the  initiation  of  an  actively  propagating  macrocrack,  (2)  fatigue  crack  propagation 
j ■ under  constant  maximum  load,  and  (3)  deviations  in  propagation  behavior  (acceleration  and/or 

I retardationl  caused  by  major  load  excursions. 

^ Kngine  hardware  operates  under  complicated  stress-time-temperature  spectra,  but  labora- 

i tory  testing  must  be  done  at  selected  conditions  because  of  cost  and  time  limitations.  To  describe 

1 crack  propagation  at  conditions  where  test  data  does  not  exist,  an  interpolative  crack  propagation 

model  is  necessary. 

This  report  presents  a demonstration  of  an  interpolative  model  for  crack  propagation 
behavior  at  elevated  temperatures,  which  is  based  on  the  hyperbolic  sine.  Propagation  life 
predictions  using  this  model  are  made  for  two  crack  geometries  tested  under  stress-temperature- 
time spectra  representative  of  advanced  gas  turbine  disk  operating  conditions. 

I 

. All  fracture  mechanics  evaluations  were  performed  on  an  advanced  nickel-base  turbine  disk 

1 alloy.  (lATORIZEI)'”  IN-KX),  used  in  the  FlOO  turhofan  engine.  Specimens  specifically 

designated  to  this  contract  were  machined  from  heat  BANQ-499,  hut  a significant  amount  of 
; crack  propagation  data  existed  for  this  alloy  prior  to  the  start  of  this  program  and  is  also  used  in 

the  analyses.  Heat  treatment  consists  of  solutionization  at  20,o0°F.  stabilization  at  lf!(K)°F  and 
18(X1°F.  and  precipitation  hardening  at  12(X)°F  and  14(X1°F.  Typical  chemical  composition  is 
ll.()7C-12.4rr-lH..')Co-:(.2.Mo-4..32Ti-4.9HAl-0.7HV-U.02B-().()HZr  - balance  nickel. 

1 

Tensile,  stress  rupture,  and  creep  test  results  for  two  forgings  from  this  heat  (499-A2A  and 
. 499-A2B)  are  given  in  Table  1. 


TABLE  1.  MECHANICAL  PROPERTIES  OF 
TWO  IN-100  PANCAKE  FORGINGS 


Tensile  Properties 


Disk  No. 

Temp 

(“F) 

Yield 

Strength 

(ksi) 

Ultimate 

Tensile 

Strength 

(ksi) 

EL 

% 

RA 

% 

499- A2A 

RT 

164.5 

232.4 

22.0 

22.2 

499- A2A 

1300 

157.2 

177.0 

14.0 

22.3 

499- A2B 

RT 

164.9 

232.0 

22.0 

21.5 

499- A2B 

1300 

156.0 

179.1 

14.7 

16.4 

Stress  Rupture  Properties 

Temp 

Stress 

Life 

EL 

RA 

Disk  No. 

(°F) 

(ksi) 

(hr) 

% 

% 

499-A2A 

1350 

95 

28.0 

10.6 

15.9 

499-A2B 

1350 

95 

18.0 

8.1 

15.6 

499-A2B 

1350 

95 

19.5 

8.0 

8.2 

499-A2B 

1350 

95 

19.3 

6.7 

12.9 

Creep  Rupture  Properties 

Total 

Temp 

Stress 

0.1% 

0.2% 

Life 

Disk  No. 

r-F; 

(ksi) 

(hr) 

(hr) 

(hr) 

499-A2A 

1300 

80 



175.5 

233.2* 

499- A2B 

1300 

80 

1U.5 

142.5 

143.2* 

‘Test  Discontinued. 


SECTION  II 
THE  MODEL 


An  interpolative  model  has  been  developed  for  the  analysis  of  elevated  temperature  fatigue 
crack  propagation  data.  This  model  is  used  to  describe  the  parametric  effects  of  four  fundamental 
influences  on  crack  propagation:  frequency  (i<),  stress  ratio  (R),  temperature  (T),  and  major  load 
excursion  effects. 

GENERAL  MODEL 

This  interpolative  model  is  based  on  the  hyperbolic  sine  equation, 

log  (da/dN)  = C,  sinh  (Cj  (log  (AK)  + C,)  ) + C4  (1) 

where  the  coefficients  are  simple  empirical  functions  of  test  frequency,  stress  ratio,  and 
temperature: 

C,  = material  constant 

C,  = f,  (R,  K,  T) 

C,  = f,  (C„  R) 

C.  = f,  (i^,  R,  T) 

A more  complete  description  of  the  model  is  presented  in  Reference  1.  The  salient  features 
are  given  here. 

It  has  been  shown  that  for  IN  TOO.  the  coefficients  in  equation  1 can  be  simple  empirical 
functions  of  cyclic  frequency,  stress  ratio,  and  temperature.  Elxperience  indicates  that,  for  a given 
material,  C,  can  be  fixed  without  adversely  affecting  model  flexibility.  For  IN-100,  C,  has  a fixed 
value  of  0.5.  The  coefficients  C,  and  C,  are  logarithmic  functions  of  cycle  duration,  \h,  (Figure 
1)  and  Cj  exhibits  linear  variation  with  log  (1-R)  (Figure  2).  The  coefficients  C,  and  C\  are  also 
linear  functions  of  temperatures  (Figure  3),  from  1000°F  to  1350°F.  C,  does  not  change  with 
temperature  because  the  correlation  lin,  is  vertical.  The  relationships  for  frequency,  stress  ratio, 
and  temperature  are  given  in  Figures  4,  5,  6,  and  7. 

The  simple  relationships  previously  discussed  de.scribe  crack  propagation  of  IN-100,  at  any 
stress  ratio  and  frequency  for  temperatures  between  1000°F  to  1350°F  in  air  conditions.  The 
computational  procedure  is  schematically  represented  in  Figure  8.  First  locate  the  coefficients  on 
the  1200°F,  R = 0.1  frequency  model,  position  1.  Second,  account  for  stress  rat  o effects  by 
moving  along  a stress  ratio  model  to  position  2.  Finally,  C,  and  C,  for  the  desired  temperature  are 
determined  using  flC/flT  from  the  temperature  model,  position  3. 

EMPIRICAL  SYNERGISTIC  MODEL 

The  majority  of  crack  growth  studies  are  performed  using  constant  amplitude  loading. 
These  tests  are  not  representative  of  gas  turbine  rotating  component  operating  conditions,  which 
include  complex  stress-time-temperature  histories.  It  is  expected  that  crack  growth  relationships 
will  be  complicated  by  mi.ssion  mix  cyclic  conditions  that  result  from  throttle  excursions,  disk- 
blade  interactions,  etc. 

A generalized  predictive  model  must  account  for  the  effects  of  load  sequence.  A mission  is 
first  analyzed  to  determine  if  principles  of  linear  superposition  of  damage  were  applicable.  If  they 
were  not,  then  synergism  must  be  considered.  Since  interaction  effects  are  sensitive  to  the  loading 
sequence,  accurate  crack  propagation  models  can  only  be  expected  when  realistic  load-time 
histories  are  simulated  (Reference  2). 
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Figure  4.  Effect  of  Frequency  on  Crack  Growth  in  IN-IOO  R = 0.1,  I200'‘F 
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Figure  5.  Effect  of  Stress  Ratio  on  Crack  Growth  in  IN-100,  10  cpm,  1200°F 


iri^C  YC  LE 


SPEC  NO 

‘.OOCiSiiO 

10005^7 

10uObHi4HF 

1 ‘lOOSSS 

lOClUbOIinF 

lOOOb^C'PF 


[T]  R - 0 

|T]  R - 0 
[T]  R - 0 


Figure  ft  Effect  of  Stress  Ratio  on  Crack  Growth  in  IN-100  20  cps,  I200'‘F 
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Figure  7.  Effect  of  Temperature  on  Crack  Growth  in  IN- 100,  10  cpm,  R ^ 0.1 
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Figure  8.  Schematic  Representation  of  the  Method  for  Determining 
SINH  Model  Coefficients  Representing  Any  Frequency,  Stress 
Ratio,  and  Temperature 


If  a realistic  load-time  mission  is  broken  down  into  small  segments,  the  segments  can  be 
represented  by  synergistic  empirical  models.  Linear  superposition  of  damage  for  each  mission 
segment  can  be  used  to  predict  crack  propagation  life.  Care  must  be  taken  when  dividing  the 
mission  to  ensure  minimum  segment  interaction  at  the  division  points.  For  example,  major  load 
excursions  affect  the  cycles  following  to  a greater  degree  than  they  are  affected  by  those  which 
proceed,  mission  division  should  be  immediately  prior  to  an  excursion  rather  than  after.  (See 
Figure  9,  .segments  4 and  7.) 

The  interpolative  hyperbolic  sine  model  is  ideal  for  representing  crack  propagation  data  of 
complex  missions  since  it  requires  only  a minimum  number  of  tests.  Both  acceleration  and 
retardation  of  the  crack  growth  rate  caused  by  overloads  can  be  accurately  represented.  The 
SINH  does  not  attempt  to  model  the  micromechanical  physical  deformation  near  the  crack  tip, 
but  empirically  describes  the  macroscopic  behavior  of  the  crack. 

This  model  is  discussed  in  detail  in  Section  III. 
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Figure  9.  Incrementation  of  Simulated  Turbine  Disk  Mission  for  I20G°F, 
Retardation  Model 
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SECTION  III 

MODEL  DEMONSTRATION 


Two  geometries,  surface  flaw  (Figure  10)  and  modified  compact  tension  (Figure  11) 
specimens,  were  tested  at  the  same  two  temperatures  to  demonstrate  the  interpolative  predictive 
ability  of  the  hyperbolic  sine  model  under  simulated  advanced  gas  turbine  disk  operating 
conditions.  Load  vs  time  for  the  mission  is  shown  in  Figure  12.  The  two  test  temperatures  were 
1000‘’F  and  1200°F. 

The  stress  intensity  (K)  solution  for  the  Modified  Compact  Tension  (MCT)  specimen 
(Reference  3)  is  shown  in  Figure  11.  The  surface  flaw  specimen  was  designed  by 
P&WA/Commercial  Products  Division  (Reference  4).  Figure  13  compares  the  analytically 
determined  surface  flaw  K-solution  and  the  solution  calculated  using  handbook  values 
(Reference  5)  for  a semicircular  surface  flaw  under  uniform  tension. 


All  testing  was  performed  on  servocontrolled  hydraulic  rigs.  MCT  specimens  were  tested  on 
an  MTS  rig  with  digital  computer  control  and  the  surface  crack  specimens  were  tested  on  a rig 
designed  and  built  by  P&W A/Florida,  using  DATATRAK*  controllers.  All  specimens  were 
heated  using  clamshell  resistance  furnaces. 

A mathematical  model  describing  crack  propagation  has  utilitarian  value  only  insofar  as  it 
can  be  used  in  life  prediction;  the  overall  accuracy  of  a model  can  be  measured  by  the  accuracy 
of  the  resulting  prediction.  To  provide  a basis  for  comparing  the  accuracy  of  various  life 
predictions,  a simple  correlative  parameter  is  used,  Npreo/N.ct,  the  quotient  of  predicted  and 
actual  cyclic  lives.  Ideally,  this  quotient  is  1.0  and  decimal  deviations  from  the  ideal  can  be 
quickly  interpreted  as  percent  error  of  the  prediction. 

A simple  cycle-by-cycle  (or  mission-by-mission)  integration  is  used  to  sum  the  incremental 
crack  advances,  da,  which  comprise  the  cyclic  life,  N,  (or  mission  life,  M). 


da/dN  = 

f (a.  AK  ) 

(2) 

dN  = 

da/f  (a,  AK  ) 

(3) 

N =f 

r 

dN  4 da/f  (a.  AK  ) 

(4) 

J 

a, 

Becau.se  this  integral  can  be  difficult  to  evaluate  directly,  computerized  numerical  methods  are 
used . 

In  equation  2,  f (AK)  is  represented  by  any  empirical  model  that  accurately  predicts  the 
average  instantaneous  crack  growth  rate  for  a calculated  stress  intensity.  Linear  cumulative 
damage  within  a mission  spectrum  is  determined  using  cycle-by-cycle  (or  unit  time-by-unit  time) 
integration  and  subsequent  summation  of  individual  components. 


N 


ininxt')tiA 


No 


da 

l(Ai^ 


E 


da 

Tuk,) 


E 

0 


da 

fO^ 


cvclic 


dwell 


sustained 

load 


(5) 


13 


Thickness  (B) 


K Calibration: 


K = Y — 
BW 


■/a 


for  a/w  = 0.3  0.7;  H/W  - H/W  = 0.6 

Y = f (a/w)  = [0.2960  - 1.855  (a/w)  + 6.557  (a/w)2 

10.17  (a/w)3  + 6.389  (a/w)^l  102 

Accuracy:  0.5% 

Net  Section  Stress; 

KW’/’ 


ONet 


f (a/w)  (w-a) 


"1  + 3(w+a)1 
(w-a)  J 


Figure  II.  Modified  Compact  Tension  Specimen 


A computer  program  has  been  developed  to  perform  crack  propagation  life  analyses  using 
(1)  linear  superposition  of  individual  components  within  a mission,  and  (2)  empirical  synergistic 
models.  The  program  predicts  the  life  of  specimens  tested  under  variable  amplitude  cyclic, 
cyclic-dwell,  and/or  dwell  conditions.  Currently,  a total  of  twenty  separate  segments  within  a 
given  mission  can  be  incrementally  integrated  to  predict  crack  propagation  life.  The  empirical 
inputs  describing  crack  propagation  within  the  mission  segments  are  the  interpolative  SINH 
model  coefficients. 

The  program  provides,  as  output,  computer  drawn  (GOULD)  plots  of  the  actual  a vs  N data 
(included  as  input)  and  the  predicted  a vs  N relationship.  This  allows  direct  comparison  of 
predicted  vs  actual  data.  Appendix  C shows  an  example  of  the  complete  program  input. 

Instructions  for  use  of  the  program  are  given  in  Appendix  A.  The  actual  program  is  listed  in 
Appendix  B. 

LIFE  ANALYSIS  AT  1000°F 

Figure  14  shows  the  test  mission  segmented  into  appropriate  regions  for  analysis  at  1(XX)°F, 
Two  assumptions  are  made  in  this  analysis.  First,  it  is  Msumed  that  sustained  load  crack 
propagation  does  not  occur  under  the  test  conditions.  This  is  based  on  the  observation  that  there 
was  no  significant  difference  between  10  cpm  and  2-minute  dwell  crack  propagation  rates  at 
1CKX)°F  (Reference  1).  Second,  it  is  assumed  that  no  synergism  occurs.  This  is  necessary  because 
there  are  no  data  yet  available  on  overload  effects  at  lOOO’F  for  IN-100.  Therefore,  the  predictions 
are  based  on  interpnilations  of  the  SINH  model  for  the  cyclic  portions,  assuming  sustained  load 
crack  propagation  to  be  insignificant. 


Figure  14.  Incrementation  of  Simulated  Turbine  Disk  Mission  for  I000°F 


17 


All  predicted  mission  lives,  and  a vs  N relationships  were  computed  and  published  prior  to 
testing,  in  PWA  FR-7906  Quarterly  Progress  Narrative  No.  6,  15  September  1976. 

The  SINH  equation  for  10  cpm,  R = 0.1  crack  growth  rate  at  1000°F  was  used  as  the  starting 
point  for  calculating  the  hyperbolic  sine  coefficients  for  each  segment  of  the  1000°F  mission.  It  is 
assumed  that  stress  ratio  effects  at  1000°F  are  similar  to  stress  ratio  effects  at  1200°F.  Since 
frequency  effects  on  crack  propagation  rates  would  be  less  at  1000°F  than  at  1200°F  because  of 
reduced  environmental  degradation  (oxidation,  Reference  1),  the  difference  between  1 cpm  and 
10  cpm  crack  growth  rates  at  1000°F  is  considered  to  be  negligible. 

The  prediction,  published  before  testing,  for  the  surface  flaw  specimen  (Figure  15>  was 
conservative.  The  accuracy  (Np„a/N,ctu.i  = 0.91)  of  the  prediction  is  encouraging,  since  only  two 
tests  had  been  conducted  at  1000°F  (10  cpm,  R = 0.1  and  2-minute  dwell,  R = 0.1)  prior  to  this 
mix  mission  test.  The  actual  a vs  N data  has  the  same  shape  as  the  predicted  a vs  N relationship 
indicating  both  an  accurate  K-solution  and  an  accurate  predictive  model. 

The  prediction,  published  prior  to  testing,  for  the  MCT  specimen  (Figure  16)  was 
anticonservative.  The  accuracy  was  good  (Np„d/N,c,u„  = 1.12),  and  the  shape  of  the  actual  data 
is  close  to  the  predicted  shape.  No  explanation  is  offered  for  the  anticonservative  behavior  of  the 
MCT  specimen. 

Tables  2 and  3 present  the  SINH  equation  coefficients,  and  other  parameters  for  the  surface 
flaw  specimen  and  the  MCT  specimen,  respectively. 

LIFE  ANALYSIS  AT  1200°F 

Figure  17  shows  the  test  mission  segmented  into  appropriate  parts  for  linear  superposition 
(no  synergistic  interactions)  life  analysis  at  1200°F.  The  interpolative  hyperbolic  sine  mode) 
(Reference  1)  provides  crack  growth  equations  for  exact  operating  conditions  during  each  segment 
of  the  mission.  Linear  superposition  of  the  damage  caused  by  each  part  of  the  mission  (no 
synergistic  interactions)  resulted  in  conservative  predictions  for  both  the  surface  flaw 
(Npred/Nactuai  = 0.57)  and  MCT  (Npfed/N,ctu«i  = 0.42)  specimens.  Figures  18  and  19  respectively 
illustrate  the  results,  and  Tables  4 and  5 give  the  SINH  equation  coefficients  and  other  data  for 
each  segment  of  the  mission. 

The  more  realistic  prediction  uses  a retardation  SINH  model  (Reference  6)  to  represent 
accurately  the  damage  caused  by  the  mission.  The  retardation  model  is  segmented  as  shown  in 
Figure  9 with  SINH  equation  coefficients  for  each  part  shown  in  Tables  6 and  7.  The  retardant 
or  accelerative  effect  of  overloads  (mission  parts  4,  7 and  10,  Figure  9)  has  been  modeled  for  the 
first  5 minutes  of  sustained  load  crack  growth  using  the  hyperbolic  sine.* 

Figure  20  illustrates  the  effects  of  25%  and  50%  repetitive  overloads  on  sustained  load  crack 
propagation  in  IN-100  at  1200‘’F.  Sustained  load  data  (no  overloads)  are  compared  with  sustained 
load  data  plus  a 25%  overload  every  2 minutes,  and  with  sustained  load  data  plus  a 50%  overload 
every  2 minutes. 


•While  the  mission  data  used  2-minute  dwell  periods,  it  has  been  shown  previously  (Reference  7)  there  is  no  significant 
difference  between  1-,  2-.  and  .'i-minute  constant  amplitude  tensile  dwell  cycles.  However,  there  is  a difference  between 
,5-minute  and  10-minute  dwells.  Therefore,  it  is  postulated  that  there  is  a critical  incubation  period  (between  5 ard  10 
minutes)  before  oxidation  becomes  the  crack  propagation  rate  controlling  mechanism. 


18 


HfiLF  CRflCK  LENGTH,  INCHES 


3F  E I M E N I'i  IJ  B E n 

651 

SiJF.'^RCE 

CFFiC^ 

0 . C h S 

1 ‘■:r  K'-S 

••IPX  LiThB 

50 . '.''00 

t*‘  I S 

( c i li"  E R H 1 Li n u 

1 ijLIij 

L‘  L L'  r i L.  L.  I' 

ETRESS  RRTin 

w T c n 

1 ' 1 - 

FHEQIJENCl 

LijhE 

Mi  55  ION 

TOTAL 


"1 1 SS 


S 


Figure  15.  Surface  Flaw  Specimen  Crack  Growth  F*rediction,  Interpolative  Hyperbolic  Sine 
Model,  No  Synergistic  Interaction,  Linear  Superposition 
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Time  Minutes 


Figure  17.  Incrementation  of  Simulated  Turbine  Disk  Mission  for  1200°F,  Linear 
Superposition 


TABLE  2.  SPECIMEN  651  SURFACE  FLAW  SPECIMEN 
1000°F 


Width  = 1.174  Thicknefts  = 

0.428 

Initial  Crack  Leniitth  = 0.030 

SI.\H  Coefficients 

Max’ 

Stre.ss 

Sfission 

Cl  C2  Cl 

C4 

Load 

Ratio 

Part 

0.500  0.320  -1.417 

4.550 

23.000 

0.1 

1 

0.5(X)  4.151  -1.278 

4.879 

,50.000 

0.5 

2 

0..500  4.151  -1.278  - 

4.679 

44. .500 

0.5 

3 

0.500  4.151  -1.278  - 

4 679 

50.000 

0.5 

4 

'All  loads  are  in  Kips. 

TABLE  3.  SPECIMEN  662  MCT  SPECIMEN  1000"F 


Width  = 2.503 

Thickneaa  * 

0.847 

Initial  Crack  Length  = 1.000 

S/.VH  Coefficients 

Max’ 

Stress 

Mission 

Cl  C2 

Cl 

C4 

Load 

Ratio 

Part 

0.500  3.320  - 

1.417  - 

4.5.50 

4.600 

0.1 

1 

0..500  4.151  - 

1.278  - 

4.679 

10.000 

0.5 

2 

0.500  4.151  - 

1.278  - 

4.679 

8.900 

0.5 

3 

0.500  4.151  - 

1.278  - 

4.679 

10.000 

0.5 

4 

'All  loads  are  in 

Kipa 
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Half  Crack  Length  - inches 


CRRCK  LENGTH,  INCHES 


SPECIMEN  NUMBER 

653 

MCT  SPECIMEN 

initihl  crock 

l.OOS 

INCHES 

MRX  LORD 

6.000 

KIPS 

TEMFERRTURE 

1200 

DEGREES  F 

STRESS  RATIO 

MIXED 

FREQUENC  i' 

MIXED 

MISSION 

A PCTIJRL  DATA 

815  MISSIONS 

PREDICTED  339  HISSIDNS 


1 2 3 

10  10"  10’ 


TOTAL  MISSIONS 


Figure  19.  Modified  Compact  Tension  Specimen  Predicted  Crack  Growth,  Interpolative 
Hyperbolic  Sine  Model,  No  Synergistic  Interaction,  Linear  Superposition 
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TABLE  4.  SPECIMEN  650  SURFACE  FLAW  SPECIMEN 
1200°F,  UNEAR  SUPERPOSITION 


Width 

= 2.509 

Thickness  = 

0,a51 

Initial  Crack  l^enKth  ~ 1.0083 

SISH  Coefficients 

Max* 

Stress 

Mission 

Cl 

C2 

C3 

C4 

Load 

Ratio 

Part 

O.WX) 

0.938  - 

■1.397  - 

4.156 

2.760 

0 20 

1 

0.500 

4.297 

1.479  - 

2.519 

2.760 

1.00 

2 

0„5(X) 

4.116 

1 .305  - 

4.241 

6.000 

0.46 

3 

0.500 

4.297 

1.479  - 

2.519 

4.680 

1.00 

4 

0..500 

4.169 

1 .277  - 

4.267 

5.340 

0..52 

5 

0.500 

4.297 

1.479  - 

2.519 

4.680 

1.00 

6 

0..500 

4.116  - 

■ 1 .305  - 

4.241 

6.000 

0.46 

7 

0..500 

4.297  - 

1,479  - 

2.519 

2.760 

1.00 

8 

•All  loads  are  in 

Kips. 

TABLE  5.  SPECIMEN  663  MCT  SPECIMEN  1200°F,  LINEAR 
SUPERPOSITION 


Width 

= 1.171 

Thickness  = 

0.449 

Initial  Crack  I.enuth  = 0.029.50 

S/.\'H  ( 

'oefficient. 

H 

Max  * 

Stress 

Mission 

Cl 

C2 

(’.3 

C4 

Load 

Ratio 

Part 

0.5(X) 

3.9.38 

-1..397 

-4.1.56 

21.344 

0.20 

1 

0,500 

4.297 

-1.479 

-2.519 

21. .344 

1.00 

2 

0.500 

4.116 

-1.305 

4.241 

46.400 

0.46 

,3 

0.5(X) 

4.297 

-1.479 

-2.519 

.36.192 

1.00 

4 

0..5(X1 

4.169 

-1.277 

4.267 

41.296 

0..52 

5 

0.500 

4.297 

- 1 .479 

-2.519 

36.192 

1.00 

6 

0..5(X) 

4.116 

1..305 

-4.241 

46.400 

0.46 

7 

0.500 

4.297 

-1  479 

--  2.519 

21. ,344 

1.00 

8 

All  loada  are  in  Kips. 


TABLE  6.  SPECIMEN  650  SURFACE 
1200°F,  RETARDATION 


FLAW  SPECIMEN 


Width 

= 1.171 

[ Thickness  = 

0.449 

Initial  Crack  I^nf^h  = 0.D295 

SISH  Coefficients 

Max* 

Stress 

Mission 

n 

C2 

C3 

C-l 

Load 

Hatio 

Part 

n.MX) 

3.958 

-1.397  - 

4.156 

21.344 

0.20 

1 

O.WX) 

4.297 

-1.479  - 

2.519 

21.344 

1.00 

2 

O.MXl 

4 116 

1 .305  - 

4.241 

46.400 

0.46 

3 

O.MX) 

4.378 

- 1 .a33 

2.998 

36.192 

1.00 

4 

O.iXX) 

4.297 

1.479  - 

2.519 

36.192 

1.00 

5 

0.500 

4.169 

-1.277 

4.267 

41.296 

0.52 

6 

0.5(X1 

4.292 

-1..565  - 

2.807 

36.192 

1.00 

7 

0.500 

4.297 

1.479 

2.519 

36.192 

1.00 

8 

O.MX) 

4.116 

- 1 .305  - 

4.241 

46.400 

0.46 

9 

0..500 

4.297 

-1.479 

2.519 

21 .344 

1.00 

10 

•All  loads  are 

in  Kips. 

TABLE  7. 

SPECIMEN 

TARDATION 

65.3  MCT 

SPECIMEN 

1200°F,  RE- 

Width 

2.509  Thickness  = ' 

0.651 

Initial 

Crack  Kenffth  = 

1 .(mi 

•S7.VH  < 

Coefficients 

.Wax- 

Stress 

Mission 

(1 

C2 

C.l 

C4 

'l.nad 

Ratio 

Part 

0..500 

3.9:w 

-1.397  - 

4.1.56 

2.760 

0.20 

1 

0..500 

4,297 

1.479  - 

2.519 

2.760 

1 00 

2 

0..500 

4 116 

- 1 ..305 

4.241 

6.1XX1 

0.46 

3 

0..5(X1 

4.378 

-1.633  - 

2.998 

4 680 

1.00 

4 

0.500 

4 297 

-1.479  - 

2.519 

4.680 

1.00 

5 

0.500 

4.169 

- 1.277  - 

4 267 

5..340 

0.52 

6 

0.500 

4. ‘292 

- 1 .565  - 

2.807 

4 680 

1.00 

7 

0.500 

4.297 

-1.479  - 

2.519 

4.680 

1.00 

8 

0.500 

4.116 

1 ,305  - 

4.241 

6.000 

0.46 

9 

0..500 

4.‘297 

-1.479  - 

2.519 

2.760 

1.00 

10 

*A1I  loads  are 

in  Kipa. 
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Fiuurr  20.  Effect  of  25'i  and  50'1  Multiple  Overloads  on  Sustained  Load  Crack  Growth  in 
IN-KX),  12(X)°F 
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The  following  equations  detail  the  relationships  between  overload  ratio  (OLR)  and  the 
resulting  SINH  coefficients  which  model  crack  growth  rate. 

The  model  used  to  determine  synergistic  SINH  equation  coefficients  for  mission  parts  4,  7 
and  10  (Figure  9)  is  defined  as  follows: 


For  1.0  < OLR  < 1.25 


For  1 .25  < OLR  < 1 .5 


C,  = 0.5 

C,  = -0.036  (OLR)  ^ 4.333 
C,  = -0.612  (OLR)  - 0.867 
C.  = -2.056  (OLR)  - 0.463 

C,  = 0.5 

Cj  = 3.016  (OLR)  + 0.518 
C,  = -0.04  (OLR)  - 1.582 
C.  = 1.25  (OLR)  - 4.598 


The  SINH  equation  coefficients  for  the  other  mission  parts  are  for  exact  conditions  using  the 
interpolative  SINH  model  developed  in  Reference  1. 

The  retardation  prediction  for  the  surface  flaw  specimen  is  very  accurate 
(N'pred/N’„,.,UH]  = 0.98)  as  shown  in  Figure  21.  The  actual  data  has  the  shape  predicted,  indicating 
an  accurate  K-solution  and  an  accurate  predictive  synergistic  model. 

The  MCT  specimen  synergistic  prediction  was  more  accurate  (Np,ed/  N,du»i  = 0.69)  than 
the  linear  superposition  prediction,  but  it  was  still  conservative  (Figure  22).  The  shape  of  the 
actual  a vs  N data  is  slightly  different  from  the  prediction.  A check  was  made  to  determine  if  this 
discrepancy  could  be  due  to  the  test  load  being  much  higher  than  precrack  load.  The  higher  load 
could  have  created  a reduction  in  crack  curvature  causing  apparent  accelerated  crack  growth  at 
the  beginning  of  the  test.  Starting  the  prediction  at  a longer  crr'^k  length,  after  the  precrack 
effects  have  disappeared  (17  missions),  matches  the  predicted  shape  to  the  actual  data  in  the 
early  part  of  the  teat  (Figure  23),  but  there  is  apparently  more  retardation  at  higher  stress 
intensities  than  predicted  by  the  model. 

I he  accuracy  of  the  M('T  K-solution  (as  well  as  the  computer  program  accuracy)  was 
checked  by  predicting  the  life  of  a simple  cyclically  loaded  specimen.  An  MCT  sjx'cimen  was 
tested  at  1 tOO°F.  R 0.3,  at  2 cpm.  and  the  actual  da/dn  vs  AK  data  was  fit  with  the  hyperbolic 
sine  etiuation.  This  regressed  e(|uation  was  then  integrated  by  the  life  analysis  program  to 
deterntine  if  the  actual  a vs  N data  could  be  reproduced,  indicating  an  accurate  K-solution. 
Figure  24  shows  the  predicted  data  matches  the  actual  data  very  well.  This  indicates  that  the 
inaccurate  mission  mix  prediction  for  the  MCT  specimen  was  not  caused  by  an  inaccurate  K- 
solution. 


Figure  25  compares  the  1(X)0°F  linear  superposition  surface  flaw  prediction  with  the  120fl°F 
linear  superposition  surface  flaw  prediction  and  the  1200°F  retardation  model  surface  flaw- 
prediction.  F.ach  prediction  used  the  same  load  and  the  same  starting  crack  lengths. 


Specimen  Numoer 
Initial  Crack 
Max  Load 
Temperature 
Stress  Ratio 
Frequency 


650  Surface  Flaw 
0.059  Inches 
46.400  KIPS 
1200  Degrees  F 
Mixed 

Luke  Mission 


Figure  21.  Surface  Flaw  Specimen  Crack  Growth  Prediction,  Interpolative  Hyperbolic  Sine  Model. 
Retardation  Model 
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CRnCK  LENGTH.  INCHES 


SPECIMEN  NUMBER 

653 

INITIAL  CRACK 

1.008 

MAX  LOAD 

6.000 

TEMPERATURE 

1200 

STRESS  RATIO 

MIXED 

FREQUENCY 

LUKE 

MCT  SPECIMEN 

INCHES 

KIPS 

DEGREES  F 
MISSION 


flCTUHL  DflTR  815  MISSION 
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Figure  22.  Modified  Compact  Tension  Specimen  Predicted  Crack  Growth,  Interpolative 
Hyperbolic  Sine  Model,  Retardation  Model 
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INITIAL  CRflCh 
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TEHPEP.HTLIRE 
STRESS  RATIO 
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1 .0D4 
6.000 
1200 
M 1 ED 
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MISSION 
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PREDICTED 
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TIAL CRACK  TO  AVOID  PRECRACK 
EFFECTS,  NOTE  THAT  THE  PREDICTED 
SHAPE  STILL  DOES  NOT  MATCH  THE 
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LENGTHS,  (i.e,,  THE  MODEL  DOES  NOT 
PREDICT  ENOUGH  RETARDATION). 
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Figure  23.  Modified  Compact  Tension  Specimen  .r^edicted  Crack  Growth,  Interpolative 
Hyperbolic  Sine  Model,  Retardation  Model 
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CRACK  LENGTH,  INCHES 


E 

[ 


SPECIMEN  NUMBER  696  MCT  SPECIMEN 

INinr.L  CRfiCK  0.38U  INCHES 

MRX  LORD  1. 195  KIPS 

TEMPERRTURE  1100  DEGREES  F 

STRESS  RATIO  0.3 

FREQUENCT  2 CPM  MISSION 


A ACTUAL  DATA  1H510  MISSIONS 
PREDICTED  19613  MISSIONS 


1 ^ 3 L 

10  10-  10  10 

TOTAL  MISSIONS 


Figure  24.  Modified  Compact  Tension  Specimen  Predicted  Crack  Growth,  Interpolative 
Hyperbolic  Sine  Model,  No  Synergistic  Interaction,  Linear  Superposition 


HALF  CP.HLK  1 ENGTH,  TNLIIL* 


Prediction  Description 

1 1200“F,  Linear  Superposition, 

No  Synergistic  Interaction 

2 1200°F,  SInh  Retardation  Model 

3 1000°F,  Linear  Superposition, 

No  Synergistic  Interaction 


Figure  2ly.  Comparison  of  1000° F and  I200°F  Linear  Superposition  Predictions  With  the  1200°F 
Retardation  Prediction 
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SECTION  IV 
CRITIQUE 


The  hyperbolic  sine  model  has  been  demonstrated  to  provide  accurate  life  predictions  for 
crack  propagation  at  elevated  temperatures  under  complex  mission  spectra.  The  accuracy  of  the 
predictions  for  the  surface  crack  specimens  (Np„<j/N,<.,uai  = 0.91  and  0.98)  was  better  than  for  the 
modified  compact  tension  specimens  (Np„d/N,rtj,i  = 1.12  and  0.69). -This  result  is  interesting 
since  the  bulk  of  the  data  used  in  model  development  came  from  MCT  specimens  (Reference  1). 

Figures  26  and  27  are  macrofractographs  of  the  fracture  surfaces  of  the  surface-crack 
specimens  tested  at  1000°F  and  1200°P',  respectively.  The  cracks  behaved  well  and  the 
assumption  of  a semicircular  crack  in  the  life  analysis  is  verified  by  the  fractographs.  As  expected, 
there  is  some  evidence  of  crack  turn  back  at  the  free  surface.  The  absolute  accuracy  of  the 
predictions  and  the  agreement  in  predicted  crack  history  (shape  of  the  a vs  N curves  relative  to 
actual  data)  leads  to  the  conclusion  that  both  the  K-solution  and  the  empirical  SINH  models  are 
accurate. 

Since  the  same  models  are  used  for  the  MCT  specimen  life  predictions,  a problem  must 
exi.st  with  test  procedures.  Figures  28  and  29  are  macrofractographs  of  the  fracture  surfaces  of  the 
MCT  specimens  tested  at  1(X)0'’F  and  1200°F.  respectively.  A shear  lip  is  observed  on  the 
specimen  tested  at  1(X)0°F.  During  constant  amplitude  load  testing  at  1(X)0°F,  0.2.')0-inch 
specimens  did  not  exhibit  shear  lips  up  to  AK  = .‘jO  KSI  The  shear  lip  observed  on  this  0.8.5- 
in.  specimen  may  therefore  be  an  unexpected  product  of  the  major  load  excursions.  Additionally, 
the  lower  temperatures  create  crack  tip  environments  that  are  more  conducive  to  mix  mode 
conditions  (Reference  6).  The  1(XX)°F  model  was  developed  with  thinner  (<  0..5  inch)  specimen 
data  and  might  be  mixed  mode  (slower  da/dN).  The  thicker  demonstration  specimen  might  be 
expected  to  exhibit  a higher  crack  propagation  rate  because  of  increased  plane  strain  constraint. 

The  specimen  tested  at  12(X)°F  has  considerable  crack  front  curvature  (lOTr).  Crack 
curvature  affects  the  accuracy  of  the  K-solution  and  therefore  the  predicted  crack  history.  This 
degree  of  curvature  was  not  generally  observed  on  thinner  (0.25  in.)  specimens  subjected  to 
mission  mix  cycling  (Reference  6).  The  curvature  observed  here  is  a result  of  the  curved  precrack 
which  resulted  from  room  temperature  (ductile)  precrack  cycling.  Room  temperature  precracking 
in  thinner  specimens  ( < 0.5  in.),  however,  produces  negligible  crack  front  curvature,  so  the  effect 
here  was  unexpected. 

The  thicker  MCT  specimens  were  machined  differently  (removed  the  chevron)  than  the 
specimens  used  to  characterize  the  material.  This  modification  was  to  reduce  curvature  during 
precrack,  but  curvature  was  still  excessive. 

The  empirical  SIN’H  model  has  been  demonstrated  to  be  an  effective  vehicle  for 
interpolative  life  prediction  under  both  simple  and  complex  loading  spectra.  Its  strength  lies  in 
its  empirical  description  of  observed  material  behavior  (crack  propagation)  rather  than  any 
quasianalytical  model  of  crack  tip  deformation  and  hypothesized  subsequent  effects  on  crack 
advance. 


SECTION  V 
CONCLUSIONS 


1.  The  interpolative  hyperbolic  sine  model  accurately  predicts  crack  propagation  life  at 
elevated  temperatures  under  complex  mission  spectra. 

2.  Linear  superposition  techniques  are  applicable  in  IN-100  at  temperatures  where  environ- 
mental degradation  is  minimal  (1000°F). 

3.  Retardation  and  acceleration  in  crack  growth  of  IN-100  due  to  major  load  excursions  should 
be  accounted  for  at  temperatures  in  the  creep  regime  (1200°F). 

4.  Further  study  is  required  to  determine  if  empirical  synergistic  models  are  generally 
applicable  to  different  component  geometries  and  crack  geometries. 

5.  Lack  of  time  prevented  running  all  four  demonstration  tests  (2  at  1000°F  and  2 at  1200°F) 
at  the  same  load.  Table  8 gives  a list  of  starting  stress  intensity  values.  A hypothetical 
comparison  is  presented  in  Figure  25. 


TABLE  8.  STARTING  STRESS  IN- 
TENSITIES 


Temperature  Starting  K 


Specimen 

Type 

°F 

ksi  Inch 

6S0 

SF 

1200 

14 

651 

SF 

1000 

15 

R52 

MCT 

1000 

,50 

65.1 

MCT 

1200 

30 

I 
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APPENDIX  A 

LIFE  ANALYSIS  COMPUTER  PROGRAM  USER  INSTRUCTION 


This  program  will  predict  crack  propagation  life  for  MCT  specimens  and  surface  flaw 
specimens  with  semicircular  flaws.  To  predict  crack  propagation  life  of  another  type  specimen, 
the  stress  intensity  solution  would  have  to  replace  that  of  the  MCT  specimen  (if  initial  crack 
length  is  greater  than  0.5  in.)  or  the  surface  flaw  specimen  (if  initial  crack  length  is  less  than  0.5 
in.).  Subroutine  life  contains  the  stress  intensity  solutions  and  the  hyperbolic  sine  equation  that 
models  crack  propagation.  The  sinh  equation  is: 

log  (da/dn)  = C,  sinh  { C,  (log  (AK)  + C,)  ) -t-  C.  (A-1) 

The  input  data  cards  must  be  in  order  as  follows: 

First  Card:  Card  columns  1-72  are  for  the  title.  This  is  in  “A"  format  under  the  variable 

name  of  TITL. 

Second  Card:  Card  columns  1-8  are  for  the  frequency.  The  designation  for  frequency  (FREQ) 

should  be  centered  in  the  first  eight  card  columns  for  output  on  the  GOULD 
plot.  FREQ  is  in  “A”  format. 

Third  Card:  Card  columns  1-8  are  for  listing  the  stress  ratio  for  GOULD  plot  output.  The 

designation  for  stress  ratio  (RR)  should  be  centered  in  the  first  eight  columns. 
RR  is  in  “A"  format. 

Fourth  Card:  Card  columns  1-52  can  be  used  for  additional  title  such  as  type  of  analysis. 

Whatever  is  input  should  be  centered  in  the  52  columns  in  “A"  format, 

^'ifth  {'ard:  Beginning  with  this  card,  the  NAMELIST  parameters  are  input  in  the  format 

shown  in  the  sample  data  set  (Appendix  C),  The  input  begins  with  &INPUT, 
then  a blank  space.  A listing  of  the  parameters  and  their  values  is  followed  by 
&END  which  ends  the  NAMELIST  input.  Commas  separate  each  parameter 
value.  Table  A-1  defines  each  NAMELIST  variable. 

Sixth  Card:  The  actual  data  is  listed  next.  A maximum  of  1(K)  data  points  may  be  input. 

Card  columns  1-10  are  for  the  crack  length  (half  crack  length  for  surface  flaw 
sj)ecimens)  and  columns  11 -20  are  for  the  number  of  missions  corresponding  to 
the  crack  length  (this  cannot  he  zero).  One  data  [Kiint  must  always  be  input, 
the  initial  crack  length  and  one  mission.  Both  variables  use  “F”  format. 

Seventh  Card:  A -1.0  is  input  in  card  columns  1-10  to  separate  data  sets  and  allow  multiple 

case  analysis. 
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TABLE  A- 1 

NAMELIST  VARIABLES 


PERIOD; 


NPER: 

COEFl: 

COEF2: 

COEF3; 

COEF4: 

PMAX: 

RATIO: 

AI; 

AF: 

TEMP: 

PKLD; 

WID: 

THIK: 

LOOPS: 


Normally  PERIOD  is  equal  to  one  for  a cyclic  portion  of  a mission,  but 
must  be  less  than  one  (unite  in  hours)  for  sustained  load  portions  of  a 
mission.  For  better  accuracy  durinR  sustained  load  portions  of  a misssion, 
PERIOD  should  be  less  than  0.1  hours.  To  shorten  PERIOD,  increase 
NPER  (PERIOD  * NPER  = total  length  of  mission  part). 

Number  of  cycles  or  time  increments  in  a mission  part. 

C,  in  the  hyperbolic  sine  crack  propagation  equation  A-1. 

C,  in  the  hyperbolic  sine  equation  A-1. 

C,  in  the  hyperbolic  sine  equation  A-1. 

C,  in  the  hyperbolic  sine  equation  A-1. 

Maximum  load  for  each  part  of  the  mission. 

Stress  Ratio  for  each  part  of  the  mission. 

Initial  Crack  length. 

Final  Crack  length  at  failure  (or  at  the  end  of  stress  intensity  solution 
validity). 

The  temperature  to  be  printed  on  the  GOULD  plot. 

The  maximum  load  for  the  whole  mission.  This  is  used  only  for  GOULD 
plot  output. 

The  width  of  the  surface  flaw  specimen  gage  section,  or  the  distance  from 
the  MCT  load  centerline  to  the  end  of  the  specimen. 

Thickness  of  the  specimen. 

The  number  of  parts  into  which  the  mission  has  been  divided. 
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APPENDIX  B 
ACTUAL  PROGRAM 


RELtASt  2.C 


MAIN 


DATE  = 7631A 


16/32/20 


LIRMt.FORT WRITTEN  bY  D.L.  SIMS  AND  C.G.  ANNIS.  LATEST 

REVISION  10/2V/76.  THIS  PROGRAM  PREDICTS  THE  LIFE  OF  MIXED 
MISSUN  SPECIMENS  FOR  PHASE  III  OF  THE  AFML  FRACTURE  MECHANICS 
CONTRACT  (MCT  AND  SURFACE  FLAW  SPECIMENS  ONLY). 


example  data  Cards  arl  as  follows 
1 LC-l  to  72  IS  FOR  THE  TITLE. 

i tC-1  TO  b IS  FOR  The  FREOUENCY  FOR  GOULD  PLOT  OUTPUT  ONLY. 

3 CC-1  TO  b IS  FOR  IHE  STRESS  RATIO  FOR  GOULD  PLOT  OUTPUT  ONLY. 

A CC-1  TO  62  IS  FUR  THE  TYPE  LIFE  PREDICTION  ( SUPERPOS IT  ION f ETC). 
f>  next  comes  namelist  data  (SEVERAL  CARDS). 

6 TnE  actual  A,N  DATA  IS  LAST  (1-100  CAROS). 

7 -l.C  MUST  BE  PLACED  BETWEEN  DIFFERENT  CASES. 

DIMENSION  A(100)<  XN(IOC),  CA(100)t  CN(100)f  COEFI(20)t  C0EF2(20) 

1 C0EF3(2C).  COEFA(20)t  TITLdB),  NPER(20),  PERIOD)  20)  ,CARO(  20  ) , 

2 PMAX(20),  FREu(2)f  RATIO(20).  RR(2).  TYPE ( 13 ) 

DATA  A,  XN.  CA,  CN/  A00*0./ 

DATA  COEFl,  C0EF2,  COE  F3  t COEFA,  PERIOD,  PMAX,  NPER/ 12C«0 . , 20  *0/ 
DATA  RATlU/20*0./ 

NAMELI  ST  /INPUT/CGEFl  ,C0£  F2,COEF3  ,COE  FA  , A I ,A  F , TEMP  ,PM  AX, 

1 WID,THIR,RA  T1  0,NPER  , PERIOD,  PKLD,  LOOPS 


C 

C 

L 

L 

C 

L 

L 

C 

C 

c 

c 

L 

L 

c 

c 


r1U=wIuTH,  THIK=(H1CKN 
title  CARO,  XN=AlTUAL 
A1  = 1NITIAL  crack  lengt 
length,  PMAX=MAX  LOAD, 
length  stored,  ta=calc 
SIDkEO,  tn=calculated 
load,  fkeq-freguency, 
FOR  g:'uld  printout  onl 

EAR! , IJOPS  = NUMBER  0 
CYCLES  UR  TIME  INCREME 
OF  LIFE  prediction  (SO 


ESS,  KATIO=STRESS  RATIO,  TITL=CCl-72  OF 
LIFE  ,.  COEF1-COEF2-COEF3-COEFA  = SINH  COEF  , 

H,  af’^critical  crack  length,  da=oelta  crack 
A = AC,TUAL  LRAC)}.  LENGTH,  CA=CALCUL ATAED  CRAC 
ulataeo  crack  length,  cn=calculated  life 
life,  dkone=k  max,  dk=oelta  k,  pkld=peak 
MISN  = NUMBER  OF  MISSIONS,  RR  = STRESS  RATI 
Y,  PERIOD  = TIME  OR  CYCLES  IN  A MISSION 
F PARTS  IN  The  mission,  NPER  = NUMBER  OF 
NTS  IN  A MISSION  PART.  TYPE  = TYPE 
PERPOSITION,  ETC). 


READ  AND  WRITE  THE  COMPLETE  DATA  SET. 
WR ITE( G,S9) 

6v  Format  ( ImI  ) 

SI  KIAD(S,S3,ENU=S5)CARD 
S3  FORMAT  (20Ar) 

WR 1 IE ( t,S7)CARD 
S7  FORMAT  ( 1X,2CA<,) 

GU  ID  SI 
SS  REWIND  S 
WRITE)  6,SV) 

10  CONTINUE 


OOCOOCIC 

00000020 

00000030 

000000*,0 

00000050 

00000060 

00000070 

ooooooeo 

00000090 

00000100 

OOOOOllO 

00000120 

OC000130 

OOOOOlsO 

00000150 

00000160 

00000170 

oooooieo 

00000190 

00000200 

00000210 

00000220 

,00000230 

000002MO 

00000250 

00000260 

00000270 

C00002BO 

00000290 

00000300 

00000310 

00000320 

00000330 

00000390 

00000350 

00000360 

K00000370 

00000380 

00000390 

000000900 

00000910 

00000920 

00000930 

00000990 

00000950 

00000960 

00000970 

00000980 

00000990 

00000500 

00000510 

00000520 

C0CC053C 

00000590 

C0C00S5C 

0OC0056C 

00000570 

00000580 
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r r r r>  o 


KELtASt  ?.0 


DATE  = 7631<» 


16/32/  20 


Initialize, 

UA  = 0. 

TN=0. 

NPTS  « C 

mount  - 0 

LESTA  = 0 
INT  = 0 
INCT  = 1 


hEAU  AND  WRITE  THE  INDIVIDUAL  PARAMETERS. 

REA015,1,ENDx9V9ITITL 
1 FORMAT (16AAI 
WR  ITH 

WRlTt(6,lVA>  TITL 
I9A  format (//IXfl6AA/» 

KEAO  The  FkEOUENCY  AND  STRESS  RATIO.  FOR  OUTPUT  ON  THE  GOULD  ONLY 
RE  A0t5  ,-,,END  = 9<»9IFRE0,  RR,  TYPE 
H FORMAT  «2A'./2A9/13A4) 

READ  The  namelist  parameters. 

RE  AOES , INPUT .ENUxVVV) 

READ  AND  WRITE  THE  ACTUAL  A VS  N DATA. 

WRIlEl  6,3V) 

iV  FORMAT E//1X, •ACTUAL  A VS  N DATA*) 

00  9 J = I ,1CC 
REAOlS.vi)  AIJ>,  XNIJ) 

V3  FORMAT (2F10.6) 

WRITE(6,V2)  AIJ),  XN(J) 

V2  FORMAT  11X,F10.7,VX,F1:  .1  ) 

IF  1 A(J  ) .LE.  0.  ) GO  TO  20 
NPTS  = NPTS  ♦ 1 
V CONTINUE 

WPIIE  SOME  OF  The  namelist  PARAMETERS. 

2C  WRITE!  6,SV) 

wR ITE( 6,2C3) WIO, THIR,A 1 

203  FORMAT  t/ZIX, 'WlOlM  = *,  FE . 6 , 6X , • T HICKNE SS  =»,F8.5,6X, 
i'lNITlAL  CRACK  LENGTH  =',F8.S) 

write  The  SINH  coefficients,  max  load  and  stress  ratio  FOR  EACH 

MISSION  PART. 

WK  ITEI  6,761  ) 

761  format 1//1AX, »SINH  COE FF ICIENTS •,16X, "MAX *,6X , 'STRESS*, 
laX  , *M1  SSION*  ) 

WRITE!  6,  762) 

762  FORMAT  !/6X  ,'C  1 • , 8X  , 'C2  • , 8X  , *C3  • , 8X  , *C *•  • , bX , *LOAD  • , 6X  , *R  ATI  O'  , 6X  , 

1 'PART'  ) 

00  b J=1 , LOOPS 

«.V2  WRITE!  6,2CO)COEF1  ! J)  , C0EF2IJ),  C0EF3!J),  C0FF4!J),  PMAXIJ), 

1 RA1I0!J),  J 

20C  Format  !//lX,<,FlC.t>,3X,F8.S,3X,Fe.3,6X,13) 
t CONTINUE 

wR1TE!6,26V)  loops 

.6V  FORMAT  !//IX, 'THE  MISSION  IS  BROKEN  DOWN  I NTO*  , I S , 3X  , ' P ART  S ' » 


00000690 
00000600 
00000610 
00000620 
00000630 
00000640 
00000650 
00000660 
00000670 
00000680 
000006  90 
00000700 
00000710 
00000720 
00000730 
00000740 
00000750 
,00000760 
00000770 
00000780 
00000790 
00000800 
00000810 
00000620 
00000830 
00000840 
00C00850 
00000660 
00000870 
C0000880 
00000690 
000009CC 
00000910 
00000920 
00000930 
00000940 
00000950 
00000954 
00000960 
00000970 
00000980 
00000990 
00001000 
00001010 
00001020 
00001030 
0000104C 
00001050 
00001060 
00001070 
ooooiceo 
00001090 
OOOOllCO 
OOOOlllC 
00001120 
00001 13C 
00001140 
00001150 


orr  o r rr  o rnor 


KtLtASl  2.C 


MAIN 


DATE  = 7631A 


16/32/20 


C initialize.  00001160 

TA=AI  , 0CC01170 

NISN=0  00001160 

0000119C 

CALCULATE  THE  LIFE  -CALL  SUBROUTINE  LIFE.  00001200 

SUbKUUTINE  LIFE  DOES  THE  INTEGRATION  CYCLE  BY  CYCLE.  00001210 

00001220 

ICl  NPAKT  = 1 00001230 

DU  155  1=1 f LOOPS  000012A0 

*•**  FIRST  PAKT  OF  MISSION  - 1ST  CALL  TO  LIFE.  **•*•*••*••••♦•••••*  00001250 

CALL  L IFE(COEFHNPART)  . C0EF2 ( NP ART  I , COEF 3 (NPART ) , C0EF61 NPART ) , 00001260 

1 AF,  lA,  TN,  WIO,  TMIK,  RATIO(NPART» , PMAXINPARTI,  CAt  CN.  00001270 

2 NPER( NPART I .NPART, PER  100 (NPART ) .MOUNT , MI SN .LEST A, INT , INCT, TI TL ,00001280 

3 AD  00001290 

00001300 

CHECK.  IF  CRACK  IS  PAST  THE  CRITICAL  LENGTH.  00001310 

IFlTA  .GE.  AF)  GO  TO  19  00001320 

00001330 

155  CONTINUE  000013*0 

00001350 

IF(TA  .LT.  AF)  GO  TO  101  00001360 

OC001370 

•••  00001360 


CALCULATE  PREDICTED  LIFE  OVER  ACTUAL  LIFE. 

1 V C = MI SN/XN(NPTS) 

C 

C WRITE  OUT  THE  RESULTS. 

WRITE)6,*13)  TITL 
*13  FOKMATdHl.lX.lbA*) 

WRITE(6,10B)  TA,  TN 

108  FORMAT ( lx, 'FINAL  CRACK  LENGTH  = • , F8. 5 , 5X , *C YC LES  TO  FAILURE  =• 
l.F  10.1  ,//) 

WRITE(6,21C)  XN(NPTS).  MISN,  U 

210  format  (//IX, 'ACTUAL  LIFE  = ' , F 7. 0, 1 OX  , 'PRE  DIC  TED  LIFE  = ', 

1 110//1  X, 'PREDICTED/ACTUAL  = ',F8.5//) 

C 

C SUBROUTINE  CMPAR  PLOTS  THE  PREDICTED  VS  ACTUAL  DATA. 

CALL  CMPAR(CA,  CN,  A,  XN , NPTS,  MOUNT,  T1TL(3),  AI  , PKLD,  TEMP, 
1 RAIIO,  FkEU,  MISN,  )OHNPTS),  RR,  TYRE) 

C 

GO  TO  10 
999  STOP 
ENU 


00001390 

00001*00 

00001*10 

00001*20 

00001*30 

00001**0 

00001*50 

00001*60 

00001*7C 

00001*60 

00001*90 

00001500 

00001510 

00001520 

00001530 

000015*0 

00001550 

00001560 

C0001570 

00001560 


o o 


RtLEASfc  i.C 


LIEfc 


□ATE 


76314 


16/32/20 


SUrtKOUTlNE  LIFE(CUEFlt  C0EF2i  C0EF3i  C0EF4,  AF.TA,  TN,  WIO, 

1 THlK,  KAIIO,  PMAX,  CA , CN,  NPER , NPART,  PERIOD,  MOUNT,  MISN, 

2 LESTA,  INT,  INCT,  TITL,  AI) 

• CALCULATE  LIFE 

DATA  01  ,02  ,03,04, DE>/  2 4. 6 , -18b  . t> ,6«  .7  , -1  01  7.  , 636  .V/ 

DATA  E 1,E2,E3,E4,E*./  7 .S  , 1 . 5 , 2 . b,  3. 5 ,4  . b/ 

DIMENSION  CAdOOl,  CN(  100  I , T1TH16) 

L DO  IHE  CYCLE  6Y  CYCLE  INTEGRATION. 

1C2  DO  169  1=1, NPER 
C 

C CHOOSE  IF  SURFACE  FLAW  OR  MCT  SPECIMEN. 

IFIAI  .LT.  0.b»  GO  TO  1333 
C 

C K-SULN  FOR  MCT  SPECIMEN. 

AUW=TA/W10 

E = 0:»A0W**E1  ♦02*AUW**E2  ♦D3*AOW**E3  ♦D‘»*A0W**E4  ♦Db*AOW*»Eb 
OR ONE = I PMAX/ » THIN* SORT (WlO) )>*F 
C 

GO  TO  1334 
C 

C R-SULN  FOR  SURFACE  FLAW  SPECIMEN. 

1333  DRUN6=  1.2063*PMAX*S0RT  ( T A ) / 1 W I0*Tm  IK  ) 

C 

1334  IFIPERIOD  .LT.  1.)  GO  TO  b 
GK  = I 1.  -«ATIU»*DKONE 

C WRITE!  6,21  AOW,F,OKONE,  DK 

. FORMAT (lX,lP4E12.bl 
GO  To  6 
b UK=UKONE 

b uAuNLG=CGLFl*SlNHtC0EF2*l  ALOG101DK)*COEF3n-*COEF4 
UAON  = 10.**OAONLG 
DA  = PERIOD  * OAON 
Ta=TA  ♦ OA 
TN=TN  ♦ 1, 

169  CL'NIINUE 
C 
C 

L COUNT  only  the  number  of  missions. 

JFINPART  .EU.  II  MISN  = MISN  ♦ 1 

c 

L 

C WRITE  OUT  ALL  The  JUNK  IN  THE  LOOP. 

IFUNCT  .LT.  I)  GO  TO  1633 
ML’UNT  = M0UNT*1 
MuOULO=MUO«MOONT,bC) 

LA (MOUNT)  = TA 
CNIMOUNTl  * MISN 
1F(MU1UL0  ,NE.  1)  GO  TO  P3b 
WKlTEl6,4l3)  TITL 

413  Format ( iHi , IX, 16A4 I 

WK1TEI6,414| 

414  FORMAT  I//1X,  •CRACK  length*, 3X, 'MISSIONS  TO  A • , 7X  , • OK  • , 7X  , 

1 'OKMAX  • ,6X1 'OA/DN*//) 

b3b  WKlTEIb,4lbl  TA,  MISN,  OK,  OKONE  , DAON 
41  S FURMAT(lX,F9.6,bX,IlC,bX,lP3El2.b) 

LESTA  = INT 


OCOOlbVO 

00001600 

CC001610 

0CCC162C 

CP0C1630 

00001640 

000016b0 

OCOC166C 

00001670 

0CC0168C 

0000169C 

00001700 

00001710 

00001720 

00001730 

00001740 

C00C17b0 

00001760 

00001770 

00001760 

00001790 

OC001600 

CC0C1810 

OOCC182C 

00001830 

C0001840 

OOOOlHbO 

C0001660 

00001870 

00001880 

00001890 

OC001900 

00001910 

00001920 

00001930 

00001940 

CC0C19bC 

00001960 

00001970 

00001980 

00001990 

C00020C0 

OC002010 

CCC02020 

00002030 

00002040 

0COO2OSO 

O0C02O60 

00002070 

00002060 

00002090 

00002100 

C000211C 

00002U0 

00002130 

00002140 

000021b0 

00002160 
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KfcLtASt  2.C 


LIFE 


DATE  = 7631A 


16/32/20 


1633  INT  = lA  * ICO. 

INCT  = INT  - LESTA 
S NFART  ‘ NPART  ♦ 1 
RETURN 
END 


00002170 

00002180 

00C0219C 

00002200 

00002210 


RELtASfc  2.0 


CMPAR 


DATE  = 7631A 


16/32/20 


SUbKUUllNE  CHPARCCA,  CN,  A,  XN,  NPTS,  MOUNT,  SPNO,  A],  PMAX, 

X TEMP,  K,  PREU,  MISN,  ALIFE,  RR , TYPE  I 

(. 

C....  SUBROUTINE  CMPAR  PLOTS  THE  ACTUAL  VS  PREDICTED  DATA. 

C 

C 

OlMENS  ION  CAI  100)  ,CN(100)  ,A(100t  ,XN(  100 ) , FREQ  ( 2 ) , RR|  2 ) , TY PE  ( 1 3) 
data  XCLNTH,NXCYCL , YCLNTM,NYCYCL/1.5,  A,  1.,  6/ 

L 

CALL  PLOTSI13.,10. ) 

L 

C DRAW  The  8.5  BY  11  RECTANGLE. 

CALL  PL01(0.,1C.,2) 

L CALL  PLUT(8.5,10.,3) 

C CALL  PL0T»B.5,0.,2 ) 

L CALL  PLOT(0.,0.,3I 

C 

C CHECKOUT  SUBROUTINE  GRAPH, 

C 

CALL  PL0T12.,0.8,-3) 

C 

C CALL  OlFFERENl  ROUTINE  IF  SPECIMEN  IS  A SURFACE  FLAW. 

IF(AI  .GT,  O.T)  GO  TO  26 

CALL  CFgRPH|CA,CN, A, XN, NPTS, MOUNT, SPNO, A I, PMAX, TEMP, R,FREO, MISN, 
1 ALIFE ,RR, TYPE ) 

GO  TO  27 
C 

it  LALL  GRAPH(CA,CN,A,XN, NPT S , MOUNT , SPNO, Al , PMAX ,TEMP ,R , FREQ ,MI SN, 

1 AL IFE ,RR,TYPE I 
27  CALL  PLOnC.  ,0.  ,99V) 

L 

RETURN 

END 


00002220 
00002230 
00002290 
000022  50 
0000226C 
00002270 
00002280 
00002290 
00002300 
00002310 
C000232C 
00002330 
00002390 
0000235C 
00002360 
00002370 
00002380 
00002390 
00002900 
00002910 
000C292C 
00002930 
00002990 
00002950 
00002960 
00002970 
00002980 
00002990 
00002500 
00002510 
00002520 
00002530 
00002590 
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KELtASt  ^.0 


AHH 


DATE 


7631  A 


16/32/70 


1 M11.N, 
..t'LOTS 


T INE 

GR 

AP 

HICA,CN  ,A,XN 

,NP 

TS 

AL  1 

FE, 

kK 

,1YPE) 

CRACK 

-L 

ENGTH, 

CYCLE 

S 1 

A 

!UN 

CAI 

1} 

,CNI  1), 

Al  11, 

XNI 

11 

lUN 

XNLUG 

1 1001  ,F 

■REU(2 

),RRI 

,N  I DATA  AND  PREDICTION. 


DATA  XCLNTh.NXCYCL  f YCLNTM.NVCYCL/l  .A,  , l.<  6/ 

.COMPUTE  THE  STARTING  INTEGER  EXPONENT  OF  TEN. 


MlNN^ALOGlCt  XNIl)  ) 
mINLN=AU)G10(  CN(1)  ) 

MINX=M1N01  MINN,  MINCN  ) 

XM  1N=M  INX 
C 

C DRAW  TrtE  LIFE  PREDICTION  SEMI-LOG  GRID. 

CALL  LPGRDIXCLNTH.NXCYCL .YCLNTM.NYCVCL ,MINX,TYPE ) 

L 

C PRINT  THE  SALIENT  PARAMETERS. 

CALL  PR INTISPNO.AI ,PMAX, TEMP.R ,FRE0,M1SN, ALIFE.RR ) 

C 

L IArE  LUGS,  COMPOTE  THE  SCALE  PARAMETERS,  PLOT  DATA  AND  PREDICTION 

C 

c Tare  logs  of  x-parameter. 

DO  1C  1=1,NPTS 
XNLOG( n = ALOG10(  XN( I ) I 
1C  CCNTINUE 


.PLOT  THE  UATA. 

XNLOGI  NPIS*n=XMlN 

xNLUG(NPTS»2(  = 1./XLINTH 

A(NPTS*l)-0.6 

AtNPT5*2l=0.2/YCLNTH 

CALL  LINE(XNLUG,A,NPTS,1,-1,2» 

. tare  logs  of  x-parameter. 

UU  2C  1=1,NCALC 

XNLOGI  1)  - ALOGlOl  CNt  1 I ) 

CUN  I INUE 

.DRAW  THE  PKEUICTIUN. 

XNLOGI NCALC*! (-XM1N 

XNL0gINCALC*2)=1./XCLNTH 

CAINCALC*! 1=0.6 

CA INCALC*2 )=C.7/YlLNTH 

LALL  L INE I XNLUG,CA ,NCA LC ,1 ,r,OI 

ke  t urn 

END 


00002SSO 
C0002S60 
0C002570 
OOC02SE<C 
0CC02SY0 
C0C02600 
P0C0261C 
CCC0262C 
00002630 
Cr0026R0 
000026S0 
C0002660 
00002670 
C00P268C 
CCCC26R0 
00002700 
CCC0271C 
O0OC272O 
00002730 
000027-.0 
.CC0">27S0 
OOOC276C 
0p00277C 
CC0C278C 
CC0027VC 
00002600 
00002810 
00002B2C 
00002830 
000028‘.0 
0COO2ESO 
CC0C266C 
00002670 
C0C026P0 
CC0C26O0 
0C002RC0 
CCOC  2<*1C 
00002020 
OOOC2N3C 
CC0C29AC 
0CC029S0 
OOCC29  6C 
00002970 
000C29FC 
CC00299C 
COCO  3000 
00003010 
00003020 
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n n 


KtLtAbt  i!.C 


LPL.RD 


DATE 


76314 


16/32/  20 


subroutine  LP(iRD(XCLNTHtNXCYCLtrCLNTM,  NYC  YCLf  minx,  type  I 
DIMENSION  TYPE113) 


PLOT  AND  annotate  A SEMI-LOG  LIEE  PREDICTION  GRID 

DATA  H ,MUVR2/C.067S,  0.044/ 
mT Im2=H*2 


DRAW  The  grid 
CALL  MT&RD2(XCLNTM,NXC 

• up.Aw  The  2-shaped  inte 

CALL  PLOTIO.,0.2,3) 
CALL  EKASE(l) 

CALL  PLOTlO.,G.b,2l 
call  e rase IOI 
CALL  PLOT(-0.3,0.3?  ,2) 
CALL  PLJT(*0. 3,0.36,21 
call  PLOI to. ,0.2,2 1 


YCL,YCLNTH,NYCYCL) 

RRUPTION. 


..WRITE  FIGURE  AND  TITLE 
1DP=6. 7 

CALL  SYMbUL(3.0,nOP»3 
1UKATURY*,0.,39,1I 
CALL  SYMbUL(3.0,T0P,KT 
I 0.,31,’l 

CALL  SYMtj0L(3.C,(T0P-T 
CALL  SYMb'JL(  3.0,  (TOP-6 
lEJICTcU  CRACK  growth', 
CALL  SYMbOL(a.O,(TuP-7 
1 ,C  . ,36 , 1 ) 

CALL  SYMb0L(3.0,( top-4 


...  annotate  The  x-axis 

YLUC--C.26 
UU  10  1=1,3 
XLOC=I *aCLNTH 
CALL  NUMdER(XLOC,YLUC, 
CALL  NUMBER (( XLOC*HT  Irt 
10  CONTINUE 


,0.,0., 


II 


....ANNOTATE  THE  Y-AXIS 

CALL  NUMbER(-C.3u,0.,M, 

DU  20  1=1,6 

YLUu  = I •YCLNTH  - H0VR2 

YVAL=0.6+I*C.2 

CALL  NUMbLR(-0.30,YLOC  ,H,YVAL,C.,  II 
ag  continue 


.LAbEL  THE  X-AXIS 
CALL  SYM3UL(3. 0000,-0. 

.LABEL  The  Y axis 
call  S YMbOL(-0.50, 1.6, 

Kt  Turn 
END 


00003030 
00003040 
00003060 
00003060 
00003070 
00003080 
000030VO 
0OCC31OC 
0CCC311O 
00003120 
0CC03130 
00003140 
00003150 
00003160 
00003170 
0000318C 
00003190 
00003200 
00003210 
00003220 
00003230 
00003240 

.•HI ,HTIM2,'PCWA  materials  AND  MECHANICS  LAB00003260 

00003260 

IM?,'AFML  contract  F33615-76-C-5097',  00003270 

00003280 

.•H|,h,'FIuURE  ',0.,9,1I  00003290 

.•HI ,H, 'MODIFIED  COMPACT  TENSION  SPECIMEN  PR000033CO 

C.,66,1)  00003310 

.•HI,H,' INTERPOLATIVE  HYPERBOLIC  SINE  MODEL '00003320 

0CCC333C 

. •HI, H, TYPE, 0., 52, II  00003340 

OCOC3350 
00003360 
00003370 
000C33P0 
00003390 
00003400 
00003410 
00003420 
OC003430 
0C0C344C 
00003460 
00003460 
00003470 
C00034EC 
00003490 
C00035CO 
00003610 
00003620 
000C3630 
00003540 
CC003660 
00003560 
00003670 
00003580 
CC00369C 


H, 10,0. ,'l , 1 I 

21 , I YLOC4Hl,H, (MINX«II,0.,-1,1I 


6, HT I M2, 'TOTAL  M1SSIUNS',0.,14,1| 
HTIM2, 'CRACK  length,  INCHE S • ,90 . , 20  I 
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r r>  r 


release  2.G 


MTl>K02 


DATE 


76  31 A 


16/32/20 


SUBROU  FINE  HTGRD2  1 XC  IN  TH  .NXCVCL  t YCLNTH  ,NrCYCL  ) 

.GENERAU  SEHl-LOG  UPEN  GRlDt  XI  LOG  > Yf  I CARTESIAN  t. 
DIMENSION  GRIUlVI  .GLNTHISil 

DATA  GRID/  0.301030,  0.477121,  0.60206&,  0.698970,  0.776131, 

1 0.643096,  0.903090,  0.934243,  1.0/ 

DATA  GLNlM/0.100,  .160,  .203,  .233,  .265,  .205,  .303, .323,  .340/ 
DAI  A H ,MUVR2/0.0B73,  0.044/ 

C 

L ElNU  X-AXISIHORIZONTAL  I CYCLE  LENGTH  AND  Y-AX  IS  ( VERT  IC  AL » MIGHT. 

YHIGHT*YCLNTM*NYCYCL 

XLGNTH=XCLNTH*NXCYCL 

C 

L DRAW  TmE  overall  X-Y  RECTANGLE,  WITH  ORIGIN  AT  (0.,  0.1. 

CALL  PLOTIC.,0.,3) 

CALL  PL0I(XLGNTH,C.,2) 

CALL  PL0TIXLGNTM,YH1GHT,2» 

CALL  PLDTIO. ,YH1GHT,2) 

CALL  PLOTIC.,0.,2) 

C 

C DKAw  TmE  HORIZONTAL  GRID  LINES. 

NLINES=NYCYCL-1 
DO  1C  I»1,NL1NES 
YLDC=1 aYCLNTH 
CALL  PLOT(0.,YLOC,3) 

C DRAW  THE  LEFT  TIC  MARK. 

CALL  PLOT IH0VR2,YLDC,2 ) 

L DRAW  The  dotted  LINE. 

Call  PLOT!  I XLGNTh-mOVR  21  ,YL0C,2,1  1 

L DRAW  THE  KlOHl  II C MARK. 

CALL  PLOIIXLGNTH.YLOC,  2) 

1C  CONTINUE 
C 

C DRAW  IHE  LOGARITHMIC  VERTICAL  LINES. 

DO  20  NC-1,NXCYCL 
NCMl  = NC-1 
DO  30  1=1,9 

XLOC  = l GRID!  II*  NCM1)*XCLNTH 
CALL  PcaTIXL0C,C.,3» 

CALL  Plot ( xloc,glnthii I ,21 
3C  continue 

L draw  the  dotted  line  at  EACH  complete  cycle. 

CALL  p lU I I XLOC, ( YH  IDHT -HUVR2I ,2,1) 

C.....L'KAw  THE  TIC  mark. 

CALL  PLOTIXLOC.YhICHT,  21 
20  CONTINUE 
C 

return 

END 


00003600 
00003610 
00003620 
00003630 
00003640 
00003650 
00003660 
00003670 
00003680 
00003690 
00003700 
00003710 
00003720 
00003730 
00003740 
C0C0375C 
00003760 
000C3770 
00003780 
00003790 
00003800 
00003810 
00003820 
C000363C 
00003840 
00003830 
00003860 
C00C387C 
00003880 
C0003890 
00003900 
00003910 
00003920 
0000 39 30 
00003940 
00003950 
00003960 
00003970 
CC0C398C 
00003990 
00004000 
00004010 
C0004020 
00004030 
00004040 
CCCC4C3C 
OCOC4C60 
coo 040 70 
00004080 
C0C04C9O 
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LtAit  2.0 


PRINT 


DATE  = 7631A 


16/32/20 


C 

C 


t 

C 


L 


5.OBR00T1NE  PR  INT  I SPND,  A1  , PMAX  , TEMP  ,R  , FREO  i MlSNt  ALIFEfRRI 

...WJULU  PRINTS  SALlENl  PARAMETERS  NEXT  TO  EACH  PLOT. 
dimension  FRECH2J  ,RR«2  » 
data  H,MTiM2/o.oe7e>,  o.i7e>o/ 

IWUA  > 2.0  * Al 

Al=2.0 

X2=3.00 

X3=3.5t 

X*=^.5C 

X">=A.  JC 

Y=7.S 

LALL  STMBOLIXl,T,H, 'SPECIMEN  NUMBER C. , 1 St  1 » 

CALL  STMBOHX2,Y,H,SPN0,0.,'.,n 

CHECK  IF  MCI  SPECIMEN  OR  SURFACE  FLAW  SPECIMEN. 

IF  t Al  .Li.  O.SI  00  TO  58 

CALL  SVMbOLIX3,Y,M,'MC1  SPECI MEN • ,0. t 12 ) 

00  TO  6b 

Sb  CALL  SYMB0LIX3, Y,H, 'SURFACE  FLAW'tO.tlZI 
6b  Y=Y-hT1M2 

CALL  SYMbOLIXltYtH, 'INITIAL  CR ACK • 1 0 . 1 1 3 1 1 I 
IFlAl  .01.  0.S1  00  TO  56 
call  NUMbEKlX2tY,M,TWOA,0.,3.n 
OU  TO  57 

56  call  NUMHERIX2tY,M,Al,C. t3,l) 

57  CALL  SYMBULlX3iY,M, 'INCHES'. 0. .6) 

Y=Y-HT IM2 

CALL  SYME0L(X1,Y,H,'MAX  LOAD* ,0. ,8 ,1 ) 

CALL  NUMBERlX2.Y,M,PMAXt0..3tll 
call  SYMB0L(X3,Y,H.'KIPS',0.,A1 
y = Y-HT  IM2 

CALL  SYMBULIXltYtH, ' TE MP ER ATORE ' , P. , U , 1 I 
call  NUMBER(X2,Y,H,TEMP,C.t-l,n 
call  SYMBUL(X3.Y,H, 'DECREES  F',0.,<J) 

Y=Y-HT 1M2 

CALL  SYMBULIXltYtH, 'STRESS  RAT lU ' , 0. , 1 2 , 1 I 
call  SYMB0UX2,  Y,H,RR,  C.  ,8,  U 
Y=Y-HT  IM2 

Call  s ymbolixi , y,h , 'fr euuency* ,c . , v, i ) 
call  SYMB0HX2,Y,M,FReu,0.,a,l> 

CALL  SYMb0LIX3,Y,H, *MI SSIUN',0. ,71 
Y=-Y-hT  IM. 

Y=Y-HT 1M2 

CALL  SYMBOL  I XI, IYyC.02  5I,H,2,C.,-I,1I 
CALL  SYMbULIX2,Y,H,*AC  lUAL  DA  T A ' , 0 . , 1 1 , 1 I 
CALL  NllMbEKIX5,Y,H,ALIFF,0.,-l,l) 

CALL  SYMLOlIXA,Y,h,'MISS10NS',0.,B) 

Y = Y-HT  IM. 

CALL  PLOT! 1.75,Y,3) 

CALL  Plot  I 2.25, Y,2 ,-I) 

CALC  SYMbOHX2,Y,H,'PRE01CTta',C.,<»,n 
CALL  NUMBER  I X5,Y,h,M1SN,C.,-1  ,1  ) 

CALL  SYMbOL(X5,Y,H,'MI  SS  I ON  S • , 0 . , 8 ) 

KEIUkN 

END 


00005100 
00005110 
OC00512P 
C 0005 130 
00005150 
00005150 
00005160 
00005170 
00005180 
0000519C 
00005200 
000052 1C 
00005220 
000052  30 
00005250 
00005250 
00005260 
00005270 
00005280 
000052 VO 
00005300 
00005310 
00005320 
00005330 
00005350 
00005350 
00005360 
00005370 
00005380 
000053VO 
00005500 
00005510 
00005520 
00005530 
00005550 
00005550 
00005560 
00005570 
00005580 
0C0055V0 
00005500 
00005510 
00005520 
00005530 
000C555C 
0CCC5550 
00CC556O 
000C557O 
C0005580 
00005590 
CCP0560C 
00005610 
00005620 
00005630 
CC005650 
00005650 
C0005660 
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KELtASt  ?.( 


LEGkW 


DATE  * 76314  16/32/20 


SUbKUUTlNt  CF6RPM(  CA.CN.A.XNtNPTStNCALCtSPNO,  AI  ,PMAX,TEMP,R,FREQt 
1 Ml  SN,  ALIFE  , RK  t TYPE  ) 

C PL0T6  LKACK-LENOTH,  CYCLES  ( AtN  I DATA  AND  PREDICTION. 

DI MENS  ION  CA( 1 » ,CN  1 1), A( n ,XN( 1 ( .TYPE  113) 

UI  MENS  1 UN  XNCOU)  IOC)  ,FRE0(2  ).RR  (2) 

DATA  XLLNTM.NXCYCL.  Y CLNTH ,NYC YCL/1 .5 . 4,  1..  6/ 

C 

C CUMPUIE  )ME  SIAkflNG  INTEGER  EXPONENT  OF  TEN. 

C 

MINN-A  LOGIC)  XN(1)  ) 

m1NlN=AL0g10)  CN)1)  ) 

M1NX=M1NC(  MINN,  MINCN  ) 

XM1N»M INX 
C 

C URA»  I ME  life  PREDICTION  SEMI-LOG  GRID. 

CALL  CFGKU)XCLNTm,NXCYCL .YCLNTH.NYCYCL .MINX.TYPE ) 

L 

C FRINl  THE  salient  PARAMETERS. 

CALL  PFINTlSPNO.AI  ,PMAX,TEMP,R,FREtl,MlSN,ALlFE,RR  ) 

L 

C IArE  lugs,  compute  the  scale  parameters,  PLOT  DATA  AND  PREDICTION 

C 

C TARt.  LUGS  OF  X-PARAMETER. 

uO  1C  1=1,NPTS 

XNlOG(I)  : ALOGIC)  XN<  1 ) ) 

1C  Continue 

c 

C PLOT  THE  DATA. 

ANLCpG)  NPTS+1  ) »XMI  N 

xNlOG) NPTS«2)»1./XLLNTH 

AINPTS ♦! )=0.0 

A(NPTS»2)=0.1/YCLNTm 

CALL  L INE ( XNLOG.A.NPTS , 1 ,-l ,2) 

C 

L lARE  LUGS  OF  X-PARAMETER. 

GO  2C  1=1,NCALC 

RNLOG)  1 ) = ALOGIC)  CN)I)  ) 

2C  LUNllNUE 
C 

C DRAM  The  PREOlCnC'N. 

XNLUG) NCALC*1 )=XMIN 

XNLUG)  NC ALC*2 ) = 1 ./XCLN  TH 

CA)NCALC«1 )=G.C 

LA )NCALL*2 )-0. 1/YLLNTH 

CALL  L INE ) XNLOG.CA ,NCA LC , 1 ,0, 0 ) 

C 

RETURN 

END 


00004670 

0C004680 

00004690 

00004700 

00004710 

00004720 

00004730 

00004740 

00004760 

00004760 

OC004770 

00004760 

00004790 

00004POO 

00004810 

00004820 

00004830 

0C'C04840 

00004860 

CCC04860 

.00004870 

0000<,8P0 

00C04890 

00004900 

00004910 

00004920 

OOOC493C 

00004940 

0C0C4950 

0OC0496C 

00004970 

0000*.980 

00004990 

OOOC6OO0 

00006010 

00006020 

00006030 

00006040 

00006060 

00006060 

00006C7C 

00006080 

00006090 

C000610C 

00005110 

00006120 

00006130 

00CC614C 
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1 


KELEASE  2.0 


CFOKl) 


DATE 


76314 


16/32/20 


SUBROOllNE  CFbROCXCLNT 
DIMENSION  1YPE(13) 


H.NXCVCL. VCLNTH.NVCYC Lt MINX t TYPE  I 00005150 

00005160 

C0C05170 

MI-LOG  LIFE  PREDICTION  GRID.  00005180 

00005190 

.044/  00005200 

00005210 

00005220 

00005230 

YCL.YCLNTH.NYCYCL)  00005240 

OC005250 

RRUPTIUN.  00005260 

00005270 

0000528C 

00005290 

00005300 

0C005310 

00005320 

00005330 

C0005340 

0C0C535C 

00005360 

.*M) tHTIM2t*PCWA  materials  AND  MECHANICS  LAb00005370 

00005380 

IM2,*AFML  CONTRACT  F336 15-75-C-5097 • , 00005390 

00005400 

.•HI ,H, 'FIGURE  •.0.,9,1)  00005410 

.•HI  ,H, 'SURFACE  FLAM  SPECIMEN  CRACK  GROWTH  P00005420 

00005430 

.•HI ,H,' INTLRPOLATIVE  HYPERBOLIC  SINE  MODEL • 00005440 

00005450 

.•Ml ,M, TYPE, 0. ,52,1  I C0005460 

00005470 

00005480 

00005490 

0C0C5500 

00005510 

H, 10,0. ,-l ,1  I 00005520 

21  , (YLOC*H),H,  (M1NX*I  l,0.,-l,ll  00005530 

00005540 

OC005550 

00005560 

,0.,C.,  II  C0CC5570 

00005580 

00005590 

000C5600 

,H,VVAL,C.,  II  00005610 

00005620 
OC 005630 
00005640 

6,mTIm2, 'TOIAL  MISSIONS', 0., 14,1)  00005650 

00005660 

00005670 

mTIm;,'malF  CRACK  length,  INCHES' ,90. ,25)  00005680 

C0C05690 

00005700 

00005710 


PLUr  AND  ANNOTATE  A SEMI-LOG  LIFE  PREDICTION  GRID. 

DA1A  H ,hOVR2/0.0875,  0.044/ 

HI 1m2=h»2 


DRAW  iHt  CRID 
CALL  MTGRD2(XCLNTH,NXC 

.DRAW  The  Z-SHAPED  INTE 
call  PL0TIC.,C.2,3) 
CALL  ERASE ( 1 1 
CALL  PL01(C.,0.5,2 I 
CALL  erase (C) 

CALL  PLUK -0. 3,0.35,2) 
call  PLUT(*C.3,C.35,2) 
CALL  PLOTIO. ,0.2,2) 


YCL,YCLNTH,NYCYCL) 
RRUPT ION. 


..wR ITE 

IuP=b. 
call  S 
UiHATuR 
call  S 

1 

lALL  S 
CALL  S 
IRE  bICT 
call  S 

l,f  ..35 
CALL  S 


FIGURE  AND  1 
7 

YMbOHi.O.n 
Y',C.,3V,1I 
YM8UL(3.0,TC 
0.,31 , 
yMbOLI 3.C  , I 1 
YMEUH3.0,  n 
ION', 0.. 45,1 
YMBOLI  3.0,  ( 1 
, 1 I 

YMEOLI 3.0,1  1 


...  ANNUIATE  The  X-AXlS 
YLCL»-C..5 
DU  10  1=1,3 
XLOL=I •XCCNTH 
CALL  NUMBEHt  XLUC.YLbC, 
LALL  NUMbtKI  I XLCC'HTIM 
1C  CDNIINUE 


....annuTau  The  y-axis 
call  NUMbLKI-0.30,0.,H 
DO  20  1=1,5 
YLCC  = I AYCLNIH  - HUVH2 
YVAL=I AC. 1 

CALL  NUMbtRI-C.3C  ,YL0C 
2C  CCM  INUt 

....label  The  x-AX I s 

CALL  SYMEIULI  i.0C0C,-C. 


.LABEL  The  Y axis 
CALL  SYMBCH-C.5C,  I.C, 

RE  TURN 
END 


H, 1C,G. , — 1 ,1  I 

2)  , (YLOC»H),H,  (M1NX*I  l,0.,-l,l) 


,0.,C.,  II 


,H,VVAL,C.,  1) 


6, hT I m2, 'TOIAL  MISSIONS',0.,14,1) 


hTIm;,'hALF  crack  length,  INCHES' ,90.  ,25) 


APPENDIX  C 
SAMPLE  INPUT  DATA 


S.HtClMtN 

LUKE 

MlXfcU 

LlNPUT 


UNO 


6b3  MCI  specimen  EOK  AEML  FRACTURE  MECrlANlCS. 

CENIER  FREOUENCY  IN  THE  FIRST  ti  SPACES  OF  THIS  LINE. 

CENTER  STRESS  RATIO  IN  TmE  FIRST  t SPACES  OF  THIS  LINE. 

RETARDATION  MODEL  CENTER  IN  S2  SPACES. 


period 

* 1 • f 

.0333,  1., 

.0330, 

0. 

0136 

, 1 . , 

.0330, 

.0283, 

I.. 

.0333, 

NPER 

= 1. 

1,  I, 

1, 

1, 

7, 

1, 

3, 

I. 

1, 

COEFI 

‘ c. bcr. 

r.SGO,  C.300, 

0.300, 

c. 

30  C , 

C.30C, 

0.300, 

0.500, 

c.soo. 

C.300, 

loef; 

* 3.93R, 

A.2v7,  A. lit.. 

A. 37b, 

A. 

2V7, 

A. 1 bV , 

A.2V2, 

A.2V7, 

A. lit, 

A.2V7, 

C0EF3 

'-I. 3V7i 

-1.A7V,-I.30S, 

-1.633, 

-1. 

A7V, 

-1.277, 

-1.363, 

-I.A7V, 

-I. 30Sf 

-1.a7A, 

CCEFR 

-“A. ISbf 

-2.31V, -A. 2a1  , 

-2. VVb, 

-2. 

31A, 

-A. 267, 

-2.6C7, 

-2.319, 

“A.  ..s  I f 

-2 .31V, 

PMAX  = 

2. 7 OL  t 

2.7bC,  6.000, 

A.bbO, 

A. 

68C  , 

3.3A0, 

A. 680 , 

A . 68  0 , 

6 . OC 0 f 

2 .760, 

RATIO 

- C.2, 

l.C,  C.A6, 

l.c. 

1 

.C, 

C.32, 

1.0, 

1.0, 

0.  A6, 

l.c. 

AI  = 1 

.0:8  jc. 

AF  = 1.73C0, 

TEMP  = 

1200., 

PKLO 

= 6.0, 

Wlo  - 

i.SCvC, 

THlK  = C.E310 

, 

LOOPS 

= 10, 

1 . CCb3 

1. 

1 .^0V7 

3. 

1.0233 

17. 

1 .C3AA 

2V. 

1.0A2V 

Al. 

1 . C/At  3 

1.  Y 

- « • 

1 . 03b2 

61  . 

i .CbCb 

71. 

1 63V 

61. 

1 .07CA 

VI. 

1 .0773 

ICl. 

1 .Ctl  1 

111. 

1 .0833 

121. 

1. :v33 

12  V. 

1 . :vv  / 

1 3b . 

1 . 1<"  26 

1a7. 

1.1133 

170. 

1.12  66 

1V3. 

1.1A63 

2 3C. 

1 . 16a0 

2 6 : . 

1.17V 7 

.VC. 

1 . 1 v<  3 

Jl  V. 

1.2123 

3a6. 

1 . ? I2b 

38  A, 

1.2  6V  2 

AVC  . 

1 .2873 

a7C  . 

1 .;07a 

3C(  . 

1 .32EI 

330. 

1 . 3 3 A.. 

36C  . 

1 ,376V 

3VC  . 

1 . 3VV2 

617. 

1 . a2v2 

63C  . 

1 .A7V3 

68  2 . 

1 . A81C 

t vv. 

1 .3132 

72  1 . 

1 .A3VA 

7a2. 

54 


l.i734 

73b 

1 . 

77C 

i.6173 

7b  I 

1.6^C1 

7VC 

1.6S75 

7V9 

1 .6b<»t> 

bCb 

1.7173 

bl3 

-l.C 
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